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We calculate the thermal conductivity k in magnetic fields near Hc2 for spin-triplet supercon- 
ducting states with line nodes vertical and horizontal relative to the Ru02-planes. The method for 
calculating the Green's functions takes into account the spatial variation of the order parameter 
and superconducting flow for the Abrikosov vortex lattice. For in-plane magnetic field we obtain 
variations of the in-plane k with two-fold symmetry as a function of rotation angle where the minima 
and maxima occur for field directions parallel and perpendicular to the heat flow. The amplitude 
of the variation decreases with increasing impurity scattering and temperature. At higher temper- 
atures the minima and maxima of the variation are interchanged. Since the results for vertical and 
horizontal line nodes are almost the same we cannot say which of the two pairing models is more 
compatible with recent measurements of k in Sr2Ru04. The observed four-fold modulation of k 
in YBa2Cu307_i is obtained for d-wave pairing by taking into account the particular shape of the 
Fermi surface and the finite temperature effect. The results for n for the f-wave pairing state with 
horizontal line nodes disagree in some respects with the measurements on UPts. 



I. INTRODUCTION 

The nature of the unconventional superconducting state in the layered ruthenate Sr2Ru04EI is currently of great 
interest. The original proposal of a spin-triplet p-wave pairing state with broken tim&jreversal symmetry and a 
constant gapa has to be modified because recent specific heatEl and NMR relaxation ratea experiments indicate the 
presence of line nodes in the superconducting gap. These results have lead to the proposal of new spin-triplet pairing 
states including f-wave pairing states with line nodes vertical and horizontal with respect to the 2D-planes.lj A powerful 
tool for probing the anisotropic gap structure is the thermal conductivity k in a magnetic field H because k is sensitive 
to the relative orientations of the heat current, the field, and the nodal directions of the gap. In fact, a modulation 
of K of 4-fold symmetry with rotation of an in-plane magnetic field reflects the angular positions of the vertical line 
nodes for d,j.2__y2 -wave pairing in YBa2Cuap7_5 (YBC0)J3 The observed modulation of the heat current in UPta 
with rotation of the in-plane magnetic fielcfl seems to be in accordance with a spin-triplet f-wave pairing state that 
has horizontal line nodes.Ei 

Recently, the thermal conductivity n in SpflE|U04 has been measured in magnetic fields up to Hc2 for field directions 
perpendicular and parallel to the ab-plane.cHi3 For rotating in-plane field the in-plane k exhibits small modulations 
of 2-fold and 4-fold symmetry with respect to the rotation angle a. Since the amplitudes-of these variations are much 
smaller than the calculated amplitudes for f-wave pairing states with vertical line nodes,Eil it has been suggested that 
the actual pairing state i|CL.Sr2Ru04 has horizontal line nodes like those proposed in Ref. |[ Most theories of thermal 



conductivity in low fiel 
an isolated vortex line. 



are based on the Doppler shift of the quasiparticle spectrum due to the circulating flow in 
Here we shall employ a quite different approach valid for fields near Hc2 in which both the 
effects of the supercurrent fiow and the scattering of the-fluasiparticles by the spatial variation of the order parameter 
in the Abrikosov vortex lattice are taken into account .EJ 

In Section II we present the theory and in Section III we discuss the results. The conclusions are contained in 
Section IV. 



II. THEORY OF GREEN'S FUNCTIONS AND THERMAL CONDUCTIVITY IN HIGH MAGNETIC 

FIELDS 

In the method of Ref. ^ the Gorkov integral equations for the normal and anomalous Green's functions G and 
F with kernels given by the Abrikosov vortex lattice function are solved by expanding all functions in Fourier series 
k with respect to the sum of the spatial positions and in Fourier integrals p with respect to the difference of the 
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spatial coordinates. In calculating the spatiaL averages it suffices to consider the k = Fourier component. The 
corresponding Green's functions are given by:Ej 

Gip,Lu) = [i3-^ + ^^/¥A2|/(p)|2(A/vsin^^)«;(z)]"' ; (1) 

F{p,Lj) = -^^/^A/(p)(A/vsin(?)u;(z)G(p,c^) ; (2) 

A = (2eH)-i/2; A^^jAirj^; (3) 

= exp(— z^)erfc(— iz) ; z = (5 + ^)A/usin0 ; (4) 

u = uj + i-f: 7 = r[7V(cj)/7Vo]"^ . (5) 

Here the applied field H is along the c-axis, 9 is the polar angle between p and H, v is the Fermi velocity in the 
ab-plane, and F is the normal-state impurity scattering rate. The spatial average A? of the absolute square of the 
Abrikosov vortex lattice function is approximately given for large k by the relationscil 

(AA/v)2 = (H,2 - H)/6/3^H ; [A/vf - (H.a/H) [6Al] . (6) 

Here, f3A — < lAj* > /(< |Ap >)'^ is the Abrikosov parameter and Aq is the BCS gap parameter. For 6^0 
the argument z in Eq. (Q) becomes very large and w{z) ~ ij^z which yields according to Eqs. ([l]) and (||) the 
ordinary Green's functions G and F for a gap A/(p) where /(p) gives the p-dependence. It should be pointed out 
that the self-energy term proportional to iw(z) in Eq. yields an imaginary part that corresponds to quasiparticle 
scattering by the vortex cores. For a field H in the ab-plane with angle a between H and the a-axis one has to make 
the following replacements in Eqs. (0) - (|^): 

sin6i— ^ [sin2 6isin2((/)-a)-f (v7v)2cos2 6l]^^^ ; (a = Z(H,a)). (7) 

Here is the azimuthal angle of p in the ab-plane and v' is the Fermi velocity parallel to the c-axis. |— ■ 

We use the general Kubo formula for the electronic thermal conductivity of strong-coupling superconductorsta and 
insert the Green's functions given in Eqs. and (|^). For simplicity we consider the 2D-case where p lies in the 
ab-plane and the polar angle Q — 7r/2. The integral over energy ^ is carried out by the method of residues (see 



Appendix C of Ref. 15) and yields the following result for k^x {x = a) for H in the ab-plane: 



K 



(8) 



XX 



8r2 7o ' ' ' Jo 2n Im^o 

[l-7r(AA/v)2|/(0)|2[sin((/.-a)]-2|^z;(zo)|2 



|1 + 2(AA/v)2|/(0)|2[sin(0 - «)]-2 [1 + i^zow{zo)] |^ 



where 



zo = 2(c^ + i7)(A/v)[sin(0 - a)]-' + ^V^(AA/v)2|/(0)|2[sin(</) - a)]-^w{zo) ; (9) 

zo = icj + ij + eo)(A/v)[sin(0 - a)]-' . (10) 

Here is the position of the pole of G in the upper half complex ^-plane which is determined by the transcendental 
equation (||) . For very low temperatures it is a good approximation to take the integrand of the (^-integral in Eq. (H) 
in the limit w ^ 0. Then the solution of Eq. (^) for zq becomes purel^^ imaginary: zq — ixq. By using the boundaries 
for the w-function at imaginary arguments one obtains the estimatellj 

xo = [27(A/v) + /?]/ sini^ - a) ; Im^o - 7 + /^(v/A) , (11) 

where 

P = 2(AA/v)2|/(0)|2 { [4(7A/v)2 + 4(AA/v)2|/(0)|2 +rjsm\<l> - a)]'/^ 

+2(7A/v) }"i ; {4/tt < r; < 2) . (12) 
Then the ratio of k^x to the normal thermal conductivity k„ is approximately given by 



K 



XX _ 1^'" d(f) ^ r(A/v) /A A / ^-2| j-/i,M-2a2i 



-^2cos^(</.) ; ; [i-(AA/v)-^i/(0)r^/3^ 

Kn Jo 27r KtA/v) + f3\ 

x\l + 2[sin(0 - a)r^ (AA/v)2|/(0)|2 _ (3[f3 + 2{jA/v)] }\~' . (13) 
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III. RESULTS FOR THE FIELD DEPENDENCE OF THE THERMAL CONDUCTIVITY FOR 

VERTICAL AND HORIZONTAL LINE NODES 



We consider first the f-wavc pairing state with vertical nodes at = ±7r/4 and ±37r/4 B 

d{p)^Az{p,+iPy)ipl-pl); |/(0)p = [cos(20)]2. (14) 

In Fig. 1 (a) and (b) we have plotted our results for Kxx/nn versus H/Hc2 for r/Ao=0.5 (0.2) and (3a — 1-2. Here we 
have used the relationships (||) connecting (AA/v)^ and (A/v)^ to the ratio H/He2 . The upper solid curves in Fig. 1 
refer to the case H || c which, according to Eq. (0), is obtained from Eq. ( [l3[) by replacing sm{(f> — a) by sm9 with 
6 = 7r/2. The lower dashed curves show the case H || a, i.e., a ~ 0. The solid curve in Fig. 1(a) increases almost 
linearly with H and agrees roughly with the data for perpendicular field shown in Fig. 1 (a) of Ref. ^ at the lowest 
temperature. The solid curve in Fig. 1(b) for F/Aq = 0.2 disagrees with the data because of its upward curvature. 
The dashed curves both exhibit an upward curvature near Hc2 which is however stronger in Fig. 1(b). Fig. 1(b) is in 
fact in better agreement with the data for parallel field shown in Fig. 1(b) of Ref. |^ due to the rapid increase of the 
thermal conductivity as H approaches Hc2. 

We next consider the angular variation of the thermal conductivity in a parallel field. In Fig. 2(a) we show an 
example of Hxx/f^n vs the angle a between H and the a-axis for (AA/v)=0.2 and F/Aq — 0.1. The angular variation is 
seen to exhibit two-fold symmetry proportional to cos(2a) with the minimum at a = where H is parallel to the heat 
current and the maximum at a = 7r/2. This behavior can be explained as an effect of the density of states which has 
a minimumJipr quasiparticles travelling parallel to the field and a maximum for quasiparticles moving perpendicular 
to the fieldJlj In fact, the expression | • • • p in the denominators of equations (||) and (|3|) corresponds to \N(4>, a)p. 
In Fig. 3(a) we show the relative amplitude Ak = 2[K{a = tt/2) - K{a = 0)]/[K(7r/2) + k(0)] for F/Ao = 0.1 (solid 
curve) and 0.5 (dashed curve) vs II/IIc2. One sees that this amplitude becomes small for II/IIc2 — > 1 and that it 
decreases with increasing impurity scattering rate F/Ag. However, even for F/Aq — 0.5, the amplitudes are still 
larger percentages of k than those shown in Fig. 2(b) of Ref. ^. 

We have also considered the case where the heat current is in the direction 0^ — tt/A corresponding to the data 
for q II [110] shown in Fig. 2(a) of Ref. ||. Then the directional factor 2cos^{(j)) in Eq. ( p^ has to be replaced by 
2 cos^((/)— 7r/4) = 1 + sin(2(/)). As an example we show in Fig. 2(a) Ks/k„ vs a for AA/v=0.2 and F/Aq = 0.1. We see 
again the two- fold symmetry with the minimum occuring at a = 7r/4, i.e., for H parallel to the heat current q, and 
the maximum at a = 37r/4. This behavior is basically in agreement with the data in Fig. 2(a) of Ref. However, at 
higher fields for AA/v=0.1 (see Fig. 2(b)), a small local maximum occurs at a = 7r/4 while a large maximum remains 
at a = 37r/4. The occurence of a local maximum and two neighboring minima is plausible since the quasiparticles 
are now travelling mainly along the node of the gap where the effect of the field on the density of states vanishes. 
This behavior seems to differ from the data in Fig. 2(a) of Ref. ^ at higher fields where a small component of 4-fold 
symmetry with minimum for H || q has been extracted (see Fig. 3 of Ref. H). The order of magnitude of the relative 
amplitudes of the calculated oscillations for the heat current along the node of the gap is about the same as that 
shown in Fig. 3(a) where the heat current is directed along the antinode of the gap. 

In Fig. 2(b) we have plotted our results for Hs/nn vs a for the case where the factor 2cos^((/i) in Eq. (|l3|) has 
been replaced with 1. The resulting curve is seen to have 4- fold symmetry in a with minima at a = and a — 7^/2 
and maxima at a = 7r/4 and a = 37r/4. This form is similar to the 4- fold symmetry of k observed-experimentally in 
YBCOJa The explanation may be the following: According to the theory of thermal conductivity,lla the conductivity 
Kjj along axis i contains the factor v? where Vj is the component of the Fermi velocity along the i-axis. For a circular 
Fermi line this yields a factor proportional to cos^(0 — However, for YBCO the Fermi line is roughy a square and 
thus the component v^^ along the b-axis is approximately constant along each side of the square. It is interesting that 
the data for Sr2Ru04 yield a small component of 4-fold synupetry which is similar to that shown in Fig. 2(b). This 
has been attributed however to the 4- fold anisotropy of lic2^ For lower fields we find an interchange of the maxima 
and minima of the 4- fold variation which is in agreement with the results in Ref. Ill (see Fig. 2(a)). 

We turn now to the f-wave pairing state with vertical nodes where the term {ptrT^ij) i^i Eq. ( |l4| ) is replaced by PxPy 
and thus the squared amplitude [cos(2(/))]^ in Eq. (|l^ ) is replaced by [sin(2(/))]^ Jj The resulting thermal conductivity 
is simply obtained from Eq. (|l3|) for state (14) by making a variable transformation (f) — (f>' — tt / A \\i the integral 



which yields the new field direction angle a' = a + tt/A and the heat current direction 0^ ~ 7r/4. Thus we obtain for 
the new state the same function of angle a' as the function of a in Fig. 2(a) (solid curve) shifted by 7r/4. Vice versa, 
for the heat current along the node of the old state the variable transformation cf) — 0' + 7r/4 yields for the new state 
the same functions of the field angle a' = a — 7r/4 as in Fig. 2(a) (dash-dot curve) shifted by — 7r/4 where the heat 
current is now along 0^ = 0. The latter form of K{a) agrees with that which has been obtained by the Doppler shift 
method for low fields in Rcf. ^ 
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The authors of Ref. ^ suggested that a superconducting state with horizontaUine nodes is a better candidate for 
explaining their thermal conductivity data. We therefore now consider the stateu 

d{p) ^ Az{px + ipy)[cos{cpz) + ao] ; |/p = [cos(cp^) + ao]^ . (15) 

where ag < 1. Now the expression for Kxx/K,n in Eq. ( p^ ) has to be integrated over p^. This leads to an additional 
integral over the new variable 6 = cpz from — tt to +tt of the expression in Eq. with the squared gap amplitude 
replaced with |/p = (cos(0) + cq)^. For H || c the terms sin(0 — a) in equations (O) and (O) have to be replaced 



by 1 which yields for oq = the same result for K^x/f^n as that for state (14). For H lying in the ab-plane one has 
to carry out the double integral over 9 and (j). In Fig. 4(a) we show Kxx/i^n vs H/Hc2 for H || a and parameter 
values r/Ao=0.5 (0.1) and (3a = 1-2. The results are seen to be very similar to the dashed curves in Figs. 1(a) and 
1(b) for the state (^4). The modulation of Hxx/^n with rotation angle a has two-fold symmetry with a minimum 
at a = and a maximum at a = tt/2 as can be seen in Fig. 4(b). The plots of Kxx/i^n vs a are very similar to 
those for state ( p^ (see Fig. 2(a), solid curve). In Fig. 3(b) we show our results for the relative amplitudes of the 
modulations for impurity scattering r/Ao=0.5 and 0.1. These are nearly the same as those for state ( p^ shown in 
Fig. 3(a). The main difference between the states with horizontal and vertical line nodes is that, for the former 
states, the two-fold modulation of k with rotation angle a of the in-plane field is always the same function relative 
to the direction of the heat current while, for the latter states, the modulation depends somewhat on the direction 
of the heat current relative to the vertical nodes (see Figs. 2(a) and 2(b)). This can be seen from Eq. ( p^ ) with the 
variable transformation (f> = (p' — (p^ for the heat current in the direction 0^ which leads to the same function of the 
new rotation angle a' = a + 0k since |/P = cos^ 9 does not depend on (j). Thus we see that our results for states with 
horizontal nodes, as well as for states with vertical nodes, agree roughly with the observed thermal conductivity data. 
This is because the minimum of k in Figs. 2(a) and 2(b) of Ref. || for heat current q || [110] and q || [100] occurs in 
both cases for H || q. 

We have also calculated the thermal conductivity for the following f-wave pairing state with horizontal line nodes 
which has been proposed for UPtsfl 

d{p)^{3V3/2)Azp,{px + ipy)^; [/[^ = (27/4) cos^ (?(1 - cos^ ^^)^ (16) 

The corresponding expression for Hbb/i^n is obtained from Eq. ( p3|) by replacing the directional terms sin(0 — a) due 
to the field by the expressions in Eq. (^ with v'=v. The directional terms |/|^ due to the squared gap amplitude are 
given now by Eq. ( |l6|) . Additionally, the heat flow directional term 2 cos^ in Eq. ([T^ ) is replaced by 3 sin^ 9 sin^ (p 
for heat flow along the b-axis. The integrations over the polar angle 9 and the azimuthal angle (p yield the variation 
of Kfib vs the rotation angle a of H in the ab-plane shown in Fig. 4(b). One see that the modulation of Kbb{<^) has 
two- fold symmetry in a with a minimum at H || b || q (q is the direction of the heat flow) and a maximum at 
H II a. These results are contrary to the two- fold variation of k observed in UPts which shows a maximum for the 
field parallel to the heat fiow.Q In Fig. 4(a) we show the dependence of Kbb on H/Hc2 for a — n/2, i.e., H || q || b. 
One sees that Hbb exhibits a strong upward curvature for small impurity scattering rate F/Ao = 0.1 while it becomes 
approximately linear in H for the larger scattering rate F/Aq — 0.5. According to the measurements of k for UPts, the 
thermal conductivity Kbb for the configuration H || q || b increases almost linearly with H while, for the configuration 
H II q II c, Kcc exhibits a steep increase as H tends to Hc2.Q Therefore we have also calculated for H || q || c and 
find surprisingly nearly the same field dependence as for Kbb^ i.e., in both cases k is linear in H or increases steeply at 
Hc2. Thus we see that both our results for the field dependence of k for H || b and H || c, and the variation of K{a) 
with the in-plane field rotation angle a are in disagreement with the measurements on UPta.Ll 

Finally we discuss the w-dependence of the integrand of Eq. (^) which determines the temperature dependence of 
the thermal conductivity. It is now necessary to solve Eq. (^ for zq as a function of uj and to integrate the resulting 
expression in Eq. (||) over (j). In Fig. 5(a) we have plotted the result (without the factor uj^sech'^{uj/2T) ) vsQ = w/A 
for the state (|lj) at field directions a = and a — 7r/2 for parameter values AA/v = 0.6 and F/Aq = 0.2. The results 
for UJ = agree with those obtained previously from Eq. (|l^ ) . Most interesting is the result that the curves for a = 
and 7r/2 in Fig. 5(a) cross each other at about ~ 0.7 which means that for higher frequencies the former minimum 
of the integrand of Kxx becomes a maximum and vice versa. This leads to a reduction of the field modulation at finite 
temperatures due to the factor u)"^ sech? {u / 2T) in the integrand of Eq. (8). In Fig. 5(b) we show our result for the 
temperature dependence of the relative amplitude An which is obtained by carrying out the w-integrals in Eq. (|^) 
with the integrands shown in Fig. 5(a). One recognizes that An decreasesj-apidly with increasing T/Tc and becomes 
small of the order of a few percent at the experimental values of T/T^Jj Furthermore note that Ak even becomes 
negative in the temperature range between T/Tc — 0.25 and 0.75 which means that the minimum and maximum are 
interchanged. 

We have also calculated the frequency dependence of the integrand of Kxx in Eq. (^) with the factor 2 cos^ (p replaced 
by 1 which yields the 4- fold variation of «(«) discussed above. For lower fields, or larger AA/v, the maximum and 
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minimum of K(a) occur, for w = 0, at field directions a = and a = 7r/4, respectively. In Fig. 6 we show the curves 
K{il) for a = and a = it/A for the same parameter values as in Fig. 5(a). Note that the curves cross each other 
for increasing frequency at a much lower frequency (fJ ~ 0.3) than in Fig. 5(a). This means that the maxima and 
minima of «(«) are interfihanged at a much lower temperature than that in Fig. 5(b). This is in agreement with the 
measurements on YBCO.eI 



IV. CONCLUSIONS 

In summary, we have calculated the in-plane thermal conductivity k, near IIc2 for different spin-triplet pairing states 
with line nodes in the superconducting gap which are directed perpendicular (vertical) and parallel (horizontal)to the 
ab-plane. For fields along the c-axis, k increases approximately linearly with increasing field, while, for fields parallel 
to the ab-plane, n exhibits an upward curvature as H approaches Hc2. For rotating in-plane magnetic field we obtain 
variations K,{a) of two- fold symmetry in the rotation angle a where the minimum and maximum occur for fields parallel 
and perpendicular to the heat flow. The relative amplitude Ak of the modulation decreases with increasing field, 
impurity scattering, and temperature. At higher temperatures Ak even changes sign indicating that the minimum 
and maximum of K{a) are interchanged. Since all these results are nearly the same for the spin-triplet pairing models 
with line nodes vertical or horizontal to the Ru03-.planes, we cannot decide which one of the proposed models is 
more compatible with the recent data for Sr2Ru04.EI'E£l The smallness of the observed variations of the in- plane k for 
rotating in-plane field may well be accounted for by finite temperature effects. The only marked difference between 
our results for these two types of states is that for horizontal nodes the variation k{q) always has the same form 
relative to the direction of the heat current while, for vertical nodes K,{a) is somewhat different for the heat current 
in the direction of an antinode or a node. For the heat current in the direction of a vertical node K,{a) shows a 2-fold 
symmetric modulation for lower fields (see the dash-dot curve in Fig. 2(a)). A similar form has been obtained in 
Ref. |l| for the state with squared amplitude sin^(20). We have shown in Section III that, for this state, the dash- 
dot curve in Fig. 2(a) is shifted by — 7r/4. For higher fields this modulation K{a) develops some structure for field 
directions near a = tt/A corresponding to the direction of the node (see the solid curve in Fig. 2(b)). Therefore we 
suggest more accurate measurements of K{a) for different directions of the heat flow which could distinguish between 
the two types of pairing states. It would also be interesting to measure K,{a) at higher temperatures to see whether 
the predicted interchange of minima and maxima actually takes place. 

No signiflcant anisotropy in the intepjplane k of Sr2Ru04 was observed for rotating in-plane held, perhaps indicating 
the existence of horizontal line nodes.EHl In fact, for the states (15) and (16) with horizontal line nodes, Kcc along the 
c-axis actually does not depend on the in-plane rotation angle a. 

The f-wave pairing state with horizontal line nodes, which has been discussed for Sr2Ru04 , has also been proposed 
for UPta.El We have calculated Kg with the help of the three-dimensional p-integral and find that, for relative orien- 
tations H II q II b as well as H || q || c, the field dependence of Kg is nearly the same: it exhibits a steep upward 
curvature near Hc2 for low impurity scattering rate F/Aq while it becomes more linear for larger values of F/Aq. 
These results disagree with measurements in UPts wherej-.^ is found to be almost linear in H for the configuration 
H II q II b while it increases steeply near Hc2 for H || q || 6.13 For rotation of H in the ab-plane we obtain a modulation 
K{a) of two- fold symmetry in the angle a with a minimum for H || b aiid a maximum for H || a. This result is in 
agreement with that obtained by the Doppler shift method for low fields.la However, these results are contrary to the 
observed two-fold variation of «(«) in UPts which has a maximum for the field parallel to the heat flow.0 Since these 
measurements were carried out at higher temperatures it may be that the interchange of maximum and minimum is 
due to the finite temperature effect. However, the disagreement between theory and experiment discussed above casts 
some doubt on the validity of this f-wave pairing model for UPta. 

Our calculation of k for the f-wave pairing state with vertical line nodes (squared gap amplitude cos^(20)) also 
applies to the d^i-yi -wave pairing state which is widely believed to be realized in the high-Tc (wprates. It is 
interesting that we obtain agreement with the observed four-fold symmetric variation K{a) in YBCOH if we neglect 
the directional factor arising from the squared component of the Fermi velocity in the direction of the heat flow. 
This approximation might be valid for the nearly square Fermi surface of YCBO because the Fermi velocity is almost 
constant on each side of the square. Our result that, for lower flelds, the rpjiima and maxima in ^(a) are interchanged 
is a problem that also occurs in the Doppler shift method at lower fields.ll2l Here again, the finite frequency behavior, 
as seen in Fig. 6, and consequently the finite temperature effect calculated at the end of Section III, are responsible 
for the interchange of minima and maxima in. the 4-fold variation K.{a) of the thermal conductivity which leads to 
agreement with the measurements on YBCOfl In general, this interchange of maxima and minima is quite sensitive 
to the strength of the magnetic field and to the temperature. 

Resolution of the remaining descrepancies between experiment and theory will require more measurements of the 
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thermal conductivity, in particular, for different directions of the Mat flow. I^hould be pointed out that we have 
made some approximations in our theory of thermal conductivityuJ near Hc2-lij First, we have neglected the small 
contributions of higher Landau levels N to the Abrikosov order parameter which occur for spin-triplet pairing.tZI 
Second, we have neglected higher order Eourier coefficients of the Green's functions with respect to wave vectors k 
of the reciprocal Abrikosov vortex latticellS which become more and more important as the applied field is decreased. 
At fields near Hd the complementary method of the Doppler effect for a single vortex becomes more appropriate, 
however, there the scattering of quasiparticles by a vortex core is neglected whereas it is taken into account in an 
average way by our method. 
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FIGURE CAPTIONS 

1. Thermal conductivity Ks/ku vs applied field II/IIc2 at T=0 for vertical line nodes (state |l4|) for H parallel to 
the c-axis (solid curves) and H parallel to the a-axis (dashed curves) for different impurity scattering parameters: (a) 
r/Ao=0.5; (b) r/Ao=0.2. The Abrikosov parameter is taken as /3y!i=1.2. 

2. ks/kh vs a = Z(H,a) for in-plane field rotation at T=0 for state ( p^ with vertical nodes and impurity scattering 
rate r/Ao=0.1. (a) Heat current q in the direction of a vertical antinode (solid curve), q in the direction of a vertical 
node (dash-dot curve), and no directional factor for q (dashed curve), all for gap parameter AA/v=0.2; (b) q in the 
direction of a node (solid curve), and no directional factor for q (dashed curve), both for AA/v=0.1. 

3. Relative amplitude Ak of the variation of k{q) for rotating in-plane fields vs II/IIc2 at T=0. (a) State ( p^ ) with 
vertical nodes and r/Ao=0.1 (solid curve) and r/Ao=0.5 (dashed curve), (b) State ( |l5| ) with horizontal line nodes 
and r/Ao=0.1 (sohd curve) and r/Ao=0.5 (dashed curve). 

4. Magnetic field dependence of ks for horizontal line nodes at T^^O. (a) ks/ku vs II/IIc2 for state (Esh with 



H II q II a (dashed curves), and for state (16) with H || q || b (solid curves), for 5 = r/Ao=0.1 and 0.5; (b) Ks/k„ 
vs a = Z(H,a) for state (|l5|) with q || a (solid curve), and for state ( p^ ) with q || b (dashed curve), for parameter 
values AA/v=0.2 and r/Ao=0.1. 

5. (a) Integrand in Eq. (^ for Ks/k„ (without the factor lo^ sech'^{u! /2T)) vs reduced frequency = w/A for state 
( |l^ ) and gap parameter AA/v=0.6, for field direction a — (solid curve) and a — t: /2 (dashed curve). The impurity 
scattering rate is r/Ao=0.2. (b) Relative amplitude A/c of the variation of K{a) for state (Q) vs T/Tc, for gap 
parameter AA/v=0.6 and impurity scattering rate r/Ao=0.2 (solid curve) and r/Ao=0.5 (dashed curve). 

6. Integrand in Eq. (0) for Ks/k„ (without the factor uj"^ sech'^{uj /2T)) and no directional factor for heat flow 
(corresponding to YCBO) vs reduced frequency O = w/A for state (14) and gap parameter AA/v=0.6, for field 
direction a = (solid curve) and a — 7r/4 (dashed curve). The impurity scattering rate is r/Ao=0.2. 
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Fig. 3a 



Fig. 3 
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